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Abstract —Segmentation is a difficult and challenging 
problem in the magnetic resonance images, and it 
considered as important in computer vision and artificial 
intelligence. Many researchers have applied various 
techniques however fuzzy c-means (FCM) based 
algorithms is more effective compared to other methods. 
In this paper, we present a novel FCM algorithm for 
weighted bias (also called intensity in-homogeneities) 
estimation and segmentation of MRI. Normally, the 
intensity inhomogeneities are attributed to imperfections 
in the radio-frequency coils or to the problems associated 
with the image acquisition. Our algorithm is formulated 
by modifying the objective function of the standard FCM 
and it has the advantage that it can be applied at an early 
stage in an automated data analysis. Further this paper 
proposes a center knowledge method in order to reduce 
the running time of proposed algorithm. The proposed 
method can deal with the intensity in-homogeneities and 
image noise effectively. We have compared our results 
with other reported methods. The results using real MRI 
data show that our method provides better results 
compared to standard FCM based algorithms and other 
modified FCM-based techniques. 
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I. INTRODUCTION 



MRI segmentation techniques [13, 16] are an essential 
technique for assisting an image-based diagnosis. 
Manual segmentation is a difficult and time consuming 
task, which makes an automated breast cancer 
segmentation [23] method desirable. The automated 
segmentation [19] of MR images into anatomical 
tissues, fluids, and structures is an interesting field in 
medical image analysis. Automated segmentation 
methods based on artificial intelligence techniques 
were proposed in (Clark et al., [8]; Fletcher Heath et 
al., [12]). Gering et al. [14] proposed a method that 
detects deviations from normal brains using a multi- 
layer Markov random field framework. In the last 
decades, fuzzy segmentation algorithms, especially the 
fuzzy c-means algorithm (FCM), have been broadly 
used in the image segmentation [24] and such a success 
mostly attributes to the introduction of fuzziness for the 
belongingness of each image pixel. Fuzzy c-means [4] 
allows for the ability to make the clustering methods 
able to retain more information from the original image 



than the crisp or hard segmentation methods [20]. But 
segmentation [14] effectiveness is refusing by bias 
field, induced by radio frequency used during MRI 
exam [2, 7]. The bias field (intensity in-homogeneity) 
is induced by the radio- frequency coil in MRI and is a 
major problem in computer-based analysis of MRI 
data. A wide variety of algorithm has been developed 
for intensity non-uniformity correction (Lai and Fang et 
al. [18]). The homomorphic filtering approach to 
remove the multiplicative effect of in-homogeneity has 
also been commonly used due to its easy and efficient 
implementation (Johnston et al.[15]; Brinkmann et al. 
[6]). In addition these approaches assume that the 
intensity corruption effects are the same for different 
patients, which is not valid in general (Lai and Fang et 
al.[18] ). Dawant et al. [10] developed a two-step app 
roach for estimation of bias field. In this approach first 
"reference points " are selected for at least one tissue 
class (they used white matter) throughout the image, 
then a thin- plate spline is "least-squared" and fitted to 
the reference point data. They suggest the coefficient of 
variations as a measure for the degree of restoration. 
S.Shen et al. [22] presented a method is called 
improved fuzzy segmentation algorithm to correct the 
intensity non-uniformity during segmentation. 

Although MRI images may appear visually uniform, 
such in-homogeneities can cause serious 
misclassifications when intensity-based segmentation 
techniques are used (Ahmed et al. [2]). Another 
approach based on the fuzzy c -means (FCM) (Bezdek 
et al.[3]; Bezdek and James et al. [5]) clustering 
technique has been used for image segmentation [21]. 
In this paper, we present a modified fuzzy c -means 
(FCM) algorithm for intensity in-homogeneities 
estimation and segmentation of brain MR images. The 
bias field can deal with the intensity in-homogeneities 
and Gaussian noise effectively. It is based on the 
traditional fuzzy c- means (FCM) clustering algorithm 
and does not consider the effect of neighborhood 
attraction [1] to correct the intensity non - uniformity 
[9] during segmentation. 

The paper is organized as follows. Section 2 presents 
the basics and proposed methods of this paper. Section 
3, discusses the segmentation results and compares the 
results with other reported techniques. Section 4 
concludes the paper. 
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II. PROPOSED METHOD 

A. Fuzzy C-Means 

Clustering methods can be considered as either hard 
(crisp) or fuzzy depending on whether a pattern data 
belongs exclusively to a single cluster or to several 
clusters with different degrees. Fuzzy c-means (FCM) 
[10] is an effective clustering algorithm for fuzzy 
clustering. Fuzzy C -means (FCM) clustering algorithm 
developed in the 1970s (Dunn [11]) and extended later 
(Bezdek [3], Bezdek et al [4]). The number of clusters 
is normally passed as an input parameter. The fuzzy c- 
means algorithm is based on minimization of the 
following objective function 



./, 



TL™i d * 



(i) 



where C is the number of cluster centers or data 
subsets; n is the number of data points; f is fuzzifier 
value (1 for hard clustering, and increasing for fuzzy 
clustering); m ik is the fuzzy membership value of pixel 

2 II l|2 

k in cluster I; d ik = \\x k — vJ is the Euclidean 

distance; X k is the k th data points; V. is centriod of 

each cluster. U is the fuzzy partition matrix and V is 
the matrix of prototypes of clusters. The above FCM 
algorithm uses iterative operation to obtain U and V 
and finally minimizes the objective function. 

B Background 

The observed MRI signal is modeled as a product of 
the true signal generated by the underlying anatomy, 
and a spatially varying factor called the gain field 

Y k =X k N k ,\/kU{\,2,...,n} (2) 

Where X k and Y k are the true and observed 

intensities at the k th voxel, respectively and N k is the 

the gain field. The application of a logarithmic 
transformation to the intensities allows the artifact to be 
modeled as an additive bias field [19]. 

y k = x k + o) k p k ykn{i,2,...,n} 

(3) 
where X k and y k are the true and observed log- 
transformed intensities at the k th voxel, respectively, 
(O k is the weight at the k th voxel and p k is the bias 

field at the k th voxel. If the gain field is known, then it 
is relatively easy to estimate the tissue class by 
applying a conventional intensity-based segmenter to 
the corrected data. 



C Centre Knowledge Algorithm 

Step 1: Let X ={x l9 x 29 ,x n }dR r be a data 

n 



set, where r-Dimension. Find s = 



.,m n , where m _ - 
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and mi,m 2 , 



*;, 



c be the number of 
ascending order. 



1/ , i=l,2,...n, j=l,2,..., r, 
cluster. Arrange mi's in 



Step 2: Rearrange the data matrix in respect of its 
relabelling mean value. (i.e) 

X y — [x\ , x' 2 , x' n J . Partitioning the data 

into c groups. First group contains first s data of X'. 
Second group contains second s data of X' 

(c-l) th group contains (c-l) th s data of X'. c th group 
contains remaining all elements. 

Step 3: Making a distance tables that show the distance 
between the elements within each group, (ie) If group 

Ik k A: 1 

X x 9 X 2 , X n J , the distance table is 
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Step 4: Select maximum distance from each distance 

table of groups. If D tJ is maximum distance of k th 

group , find the mean value M k of the elements x t and 
Xj. k th cluster center = M k . k=l,2,. . .,c 

D Proposed Novel FCM [NFCM] 

The new objective functions of FCM as shown below. 
Objective function 



j m =ti4 ad *A d *) 2 } + -\ i -i^ 
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where df k = \x k - vJ and d ik 
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(4) 
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= \\Xu +£-bu -v, . 



£□(0,1) and a,S>0. 

The objective function J m can be minimized in a 
fashion similar to the standard FCM algorithm. Taking 
the first derivatives o f J m with respect to Mik , v i and 
b k and by setting them to zero results in three 
estimator of U , V , and b . With these estimators we 
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can form an algorithm to compute the tissue class and 
bias field. 

a). Bias field estimation: Taking the derivative of J m 
with respect to b k and setting the result to zero we 



3 .* 



^^{x^s-h-v^-l)) 



^®>k 



h= h k 



(5) 



Differentiating the distance expression, we obtain: 

= 






bjL =b t 



(6) 

The zero-gradient condition for the bias- field estimator 
is expressed as: 

f \ 






£□(0,1) is the weight 
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b). Cluster prototype updating: Taking the derivative of 
J m with respect to v t and setting the result to zero, we 
have 

n 

Y d MTk{{a + S)x k + S(s-b k )) 

v*=^ 

(8) 
c). Membership evaluation: We compute this using 
Lagrange multiplier as shown below 

(9) 

Taking the derivative of L m with respect to Mik and 
setting the result to zero, we have: 



Mi k 
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Since we have: 









(11) 



where c is the number of clusters, m = 2 . 

III. RESULTS AND DISCUSSION 



In this work, we present proposed modified fuzzy 
clustering method for the segmentation of Tl- T2- 
weighted brain MRI of the same patient. The brain Tl- 
T2 weighted images corrupted by "Gaussian" noise for 
the purpose of experimental work (given in Fig. l(a-b). 
In nature, the MRI images typically do not suffer from 
"Gaussian" noise, and we add such type of noise just 
for the comparison of robustness to noises of proposed 
algorithm. In this section, we describe some 
experimental results to compare the segmentation 
performance of the following algorithms, i.e. Improved 
fuzzy segmentation algorithm [IFS][16], KFCM[6], 
and proposed Novel FCM. We test these three 
methods on brain MR image. Figs. 2-4(a-b) show the 
segmentation results of IFS, KFCM, and NFCM 
respectively. As shown in Figs. 3-4(a-b), neither IFS 
nor KFCM can separate the six classes, while NFCMS 
completely succeed in correcting and classifying the 
data as shown in Fig. 4a-b. From the images, we can 
see that both IFS and KFCM are affected by the noise 
badly, while NFCM nearly completely eliminate the 
effect of noise. 





Fig. 1. (a) Tl corrupted by Gaussian noise, (b) T2 
corrupted by Gaussian noise 





Fig. 2. (a) Tl Segmented by IFS, (b) T2 Segmented by 
IFS 
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Fig. 3. (a) Tl Segmented by KFCM, (b) T2 Segmented 
byKFCM 





Fig. 4. (a) Tl Segmented by NFCM, (b) T2 Segmented 
by NFCM 



Table I. Segmentation Accuracies 



No. of 
Clusters 


Silhouette 
Value 


Accuracy 


IFS 


0.45 


45% 


KFCM 


0.51 


51% 


Proposed 
NFCM 


0.73 


73% 



Table I gives the segmentation accuracy of the three 
algorithms on two different T1-T2 noisy images, where 
segmentation accuracy is defined using silhouette value 
in (Kannan S. R. [17]). These silhouette average values 
measures the degree of confidence in the clustering 
assignment of a particular observation, with well- 
clustered observations having values near 1 and poorly 
clustered observations having values near -1. From 
Table I, the best clustering validity 78 % was obtained 
for proposed method during the experimental work on 
brain image data. Further, it is clear from Fig. 4(a-b) 
that our proposed method completely succeeded in 
correcting and classifying the breast data and almost it 
eliminated completely the effect of noise in images. 
The other two IFS and KFCM techniques partially 
corrected the misclassified pixels given in Figs. 2-3 (a- 
b). As can be seen, the clustering results of our 
Modified Fuzzy C-Means clustering algorithms were 
superior to those obtained by using other two 
algorithms which we proposed. 



IV. CONCLUSION 

This paper presents a new novel fuzzy clustering 
algorithm with canter knowledge method for 
segmentation of brain Tl and T2 weighted images. 
This paper described experimental results on real MR 
images which corrupted with Gaussian noise to show 
the segmentation performance of proposed method. 
The segmentation results of proposed method have 
compared by existing methods. Further, the 
segmentation accuracy was obtained using Silhouette 
method and the proposed algorithm produces high 
segmentation accuracy than existed methods. The 
results reported in this paper show that the proposed 
novel objective function of fuzzy c-means is an 
effective approach to construct a robust image 
segmentation algorithm. 
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